Slave Boson Approach to The Neutron Scattering in YBa 2 Cu306+j, Superconductors 
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The evolution of the so-called "41 meV resonance" in the magnetic response of YBCO cuprates is 
studied with slave-boson theory for the i-t'-J-model. The resonance appears as a collective spin 
fluctuation in the d-wave superconducting (SC) state. It is undamped at optimal doping due to a 
threshold in the excitation energies of particle-hole pairs with relative wave vector (n, n) . When 
hole filling is reduced, the resonance moves to lower energies and broadens. Below the resonance 
energy we find a crossover to an incommensurate response in agreement with a recent experiment 
on YBa2Cu30s.6 • We show that dynamic nesting in the d-wave SC state causes this effect. 



PACS numbers: 71.10.Fd, 74.25.Ha, 74.72.Bk, 75.20.Hr 



One of the puzzling phenomena in copper-oxide super- 
conductors is the so-called "41 meV resonance" observed 
in inelastic neutron-scattering (INS) experiments on op- 
timally doped YBa2Cu306+j, (YBCO) compounds in the 
superconducting (SC) phase Jl]-|j| . The magnetic suscep- 
tibility x"(q, uj) shows a peak at the antiferromagnetic 
(AF) wave vector (tt, 7r) and an energy uj q =41 meV . The 
peak appears to be resolution limited, i.e., almost un- 
damped, and vanishes above T c . In underdoped YBCO 
the peak shifts to lower energies u> Q < 41 meV and 
develops some damping. It is also visible in the spin- 
gap regime above T c . Several theoretical approaches to 
the magnetic excitations have been proposed [5-15 . In 
most cases P-p^| the resonance at optimal doping is iden- 
tified as a collective spin fluctuation, which is stabilized 
through the suppression of quasi-particle damping in the 
SC state, and vanishes in the normal phase. The ques- 
tion arises how this mode evolves with decreasing doping. 
In particular it is unclear how the finite damping at in- 
termediate doping levels can be obtained, since the peak 
is expected to narrow and become a spin-wave (Bragg) 
peak at zero energy when the Neel state is reached. There 
is also experimental evidence |^,[l6| that the magnetic re- 
sponse in YBCO at energies somewhat below lvq is incom- 
mensurate. It is dominated by four peaks at q = (71", ir±5) 
and q = (n ± 6, tt) , a pattern previously known only for 
the La2- 2; Sr 2 ;Cu04+2 / family of compounds. In view of 
the new experimental data, we re-investigate the spin re- 
sponse in slave-boson mean-field theory |p|-|ll"|| , focusing 
on the evolution of the resonance with hole filling (dop- 
ing) and the crossover commensurate-incommensurate 
with variation of energy. 

We start from the t-t'-J-model for a single Cu02-layer. 
Results specific to the bi-layer structure of YBCO will be 
presented elsewhere. The model reads 
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Sums include nearest neighbor or next n.n. 

Cu-sites on a 2D square lattice. Physical operators c%& — 
b\fia , Si = 5 J2a.a> }L T ' J ' J ' f«?' are represented by auxil- 



iary ('slave') bosons bi carrying the charge and fermions 
fi a representing the spin a . In mean-field theory the d- 
wave superconducting phase is represented by condensed 
bosons bi — + (bi) and fermions in a d-wave BCS state 
with dispersion 

e(k) = — 2t[cos(k x ) + cos^j,)] — At cos(k x ) cos(k y ) — \i$ 

and gap A(k) = Ao[cos(fc x ) — cos(fc !/ )]/2 . Renormal- 
ized hopping parameters are t = xt + \ J(fl^fi+ x ]) , 
t' = xt' . The expectation values (. . .) in t and Ao = 
2J[{fitfi+xl) - (filfi+xt)} are computed self-consistently 
from minimization of the free energy. The hole concen- 
tration (doping) x sets the density of fermions (1 — x) — 

12Mi<rfi<r) and bosons ( b i) = V% ■ 
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FIG. 1. Particle-hole (ph) irreducible contribution to the 
susceptibility Eq. (|l]) . Single and double arrowed lines are 
normal and anomalous Green's functions of auxiliary fermions 
in the d-wave BCS state. Shaded areas denote the vertex 
function A = 1 + A J . A J accounts for all ladder diagrams not 
included in a simple random-phase approximation. 

The magnetic susceptibility is calculated in ladder ap- 
proximation, 





X(q,w)=x (q,^)/[l+aj(q)x (q,w)] 



(1) 



with J(q) = J[cos(q x ) +cos(q v )} and a = 1 . Fig. |l| shows 
X (q, ijS) , which is built of fermion bubbles with a vertex 
function to include all ladder diagrams. Boson excita- 
tions do not appear in ladder approximation. Therefore 
the calculation does not distinguish between the SC and 
the spin-gap phase (Ao ^ but (bi) = 0), and the persis- 
tence of the resonance above T c in underdoped cuprates 
appears quite naturally. 

If x(q, = 0) is computed from Eq. (0) , an instability 
to an (incommensurate) Neel state occurs for hole densi- 
ties x < x c ~ 0.2 . The high value of x c is an artifact of 
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the mean-field theory. Since the inclusion of fluctuations 
is beyond the scope of this work, we model the expected 
suppression [|l7],[l8| by setting a = 0.34 in Eq.(|l|), such 
that x c is reduced to an observed value 0.02 . It is as- 
sumed that the superconducting phase x > x c is stable 
when fluctuations are included, at least for low T . A 
renormalization of the mean-field parameters, e.g., Ao is 
ignored, in order to keep a <-* x c the only phcnomcno- 
logical input. The self-consistent Ao « 45meV at op- 
timal doping actually compares to experimental values 
p9| . We have checked numerically that the local mo- 
ment sum-rule is not sensitive to the choice of a . 
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FIG. 2. Imaginary part x" °f the susceptibility Eq.(jj|) 
for wave vector (ir, ir) . Main figure: superconducting state 
at T — > ; different peaks belong from left to right to hole 
concentrations x = 0.02 (scaled x 0.1), 0.025 (x 0.1), 0.04, 
0.06, 0.08, 0.12, 0.16. All peaks except x = 0.04,0.06,0.08 
are 5-functions, computed with a quasi-particle damping 
r = 10" 3 J. Inset: x = 0.06,0.12 in the SC state (full 
lines) and x = 0.12 in the normal state T >T C (dashed line). 
Here an experimental energy resolution (FWHM) of 5 meV is 
simulated via T = 0.0 1J . 

With this approach we are able to access the full range 
of hole densities 0.02 < x < 0.15 of underdoped cuprates. 
Numerical calculations are performed in the supercon- 
ducting state at low temperature T — » with parameters 

\ shows the 
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t = 2 J , t' = -0.45* for YBCO |20 

imaginary part x"(Qir, w ) of Eq7(|l) at the AF wave vec- 
tor Qtt = (tt, 7r) for several hole concentrations. Appar- 
ently x" is dominated by a sharp resonance. For x = 0.12 
it appears at an energy uiq = 0.51 J « 60meV, its resid- 
ual width is due to the small quasi-particle damping used 
in the numerical calculation. When doping is reduced to 
x « 0.08 — > 0.04 the resonance moves monotonously to 
lower energies and develops some damping. For further 
reduced hole filling, the peak becomes again resolution 
limited and eventually shifts to ujq — > when the AF 



transition is reached at x = x c — 0.02 . The peak van- 
ishes in the normal phase T > T c near optimal doping, 
see the inset in Fig. ^|. In heavily overdoped systems 
x > 0.2 the resonance vanishes even in the SC state. 

The qualitative agreement between these results and 
the experimental findings summarized in the beginning is 
quite satisfactory. The resonance position ujq rs 60mcV 
near optimal doping (x — 0.12) is not too far off the 
observed value 40 meV . We also calculated the AF cor- 
relation length Q from the equal-time correlation func- 
tion (S(q)S(— q)) . When doping is reduced, C increases 
monotonically and reaches the system size at x = x c = 
0.02 . Quantitatively £ is overestimated by a factor of 



2 compared to known values 1 22 2 
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FIG. 3. Imaginary part x°" °f the particle-hole bub- 
bles from Fig. [j] for wave vector (it, 7t) (continuous lines), 
and the inverse Stoner enhancement-factor K defined in the 
text (scaled x(— 10) , dashed lines). Main figure: shown are 
three hole fillings x in the SC state at T — > . An arrow 
points to the threshold energy Qo in the respective x°" j the 
corresponding K crosses zero at ujq nearby. x°" shows a van 
Hove singularity at 2Ao > 0.7J (a peak) and for x = 0.12 
also at Oo « 0.53J (a step). Inset: normal state T > T c for 
x = 0.12. 

It is at first sight surprising that the resonance occurs 
at energies ujq < 2Ao without significant damping from 
particle-hole (ph) excitations, since the d-wave SC phase 
has a finite density of states. Fig. || shows the imagi- 
nary part x "(Qtt,u) of the ph-bubbles x° = X ' + ix°" 
from Fig. [f] , together with the real part of the denomi- 
nator of Eq.flij), K(Q„,u) = [1 + aJ(Q 7r )x°'(Q ff , w)] . 
From the numerical calculation the vertex corrections 
have no effect on the outcome of Eq.(|l|) and are omit- 
ted here. x°" nas a mu S a P U P to an energy fig , which 
increases with doping. A pole occurs in Eq.(Q) since the 
corresponding K crosses zero at an energy luq < Qq in 
the gap. The result is a 8-\ike resonance for a hole fill— 
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ing near the AF transition (x = 0.025) or near optimal 
doping (a; = 0.12) . In the underdoped case (x = 0.06) 
we have ujq ~ , and the resonance is asymmetrically 
broadened. To explain this we note that apart from 
BCS coherence factors, x°"(q, w) ~ J2k — ^(^k)) 
at T -> 0, where Q(q ,k) = E(k) + E(k + q) and 
E(k) = v /e 2 (k) + A 2 (k). Here n(q,k) denotes the ph- 
excitation energies of fermions with relative wave vec- 
tor q . For q = this has a minimum value f2o > 
which determines the threshold in x°" • It is given by 
fio = 2\/if\Z , where /^/ is the fermion chemical potential 
and Z = 1 for k = A 2 /(4^/?) > 2 and Z = y/n- n 2 /^ 
for < k < 2 . Note that Oo , and therefore the resonance 
energy ojo do not follow 2Ao . 

The reason for the absence of damping at optimal dop- 
ing (x — 0.12) is identified in the step-like van Hove 
singularity (v.H.s.) at ui = f^o in X°" ( see Fig- §): in- 
duced by a locally flat ph-dispersion. By virtue of the 
Kramers-Kroenig transformation the real part x°' de- 
velops a sharp structure around Qq , shifting the posi- 
tion wo of the resonance (i.e., the zero crossing of K) 
well into the gap. When x is reduced, the v.H.s. and 
with it the peak structure in x weaken, and loq moves 
close to and may cross the threshold VLq to damping ph- 
excitations. At very low doping x ~> x c this trend is over- 
compensated by the monotonous increase of \° with re- 
duced x , which shifts ujq back into the gap. The increase 
of x° comes from the shrinking of the upper cutoff 2W for 
ph-excitations, with the bandwidth W « 8t ~ 8(xt+ J/8) 
of Gutzwiller renormalized fermions. It should be noted 
that the step-like v.H.s. in x°" near optimal doping de- 
pends on the coexistence of a finite Ao and effective n.n.n. 
hopping t' , which lead to a sufficiently flat f2(Q„.,k). 
Setting t' = eliminates the step, and the resonance 
from Eq.(Q) is severely broadened. In the normal phase 
Ao = no resonance appears at all. The ph-dispersion 
n(Q 7r ,k) = |e(k) | + |e(k + Q-n-) | then has a zero minimum 
value without v.H.s., resulting in a gapless and structure- 
less x°" and K (see the inset of Fig. ||). 

The approach presented here supports the understand- 
ing of the "41 meV resonance" as a collective spin fluc- 
tuation. We find that the sharp resonance is entirely 
caused by the pole in the random-phase approximation 
(RPA), which sums up spin-singlet particle-hole (ph) ex- 
citations, i.e., transversal spin fluctuations. The vertex 
function A = 1 + A' 7 in x° (see Fig. |j]) has numerically 
no effect (A J 0) in the energy and doping range con- 
sidered here. That is, we can neglect in particular the 
resonant contribution from spin-triplet particle-particle 
(pp) pairs, which are involved through the mixing with 

ph-excitations /fe + Q t/fej fk+Q w -tf-kt m the SC state 
Jl5| . Our view on the neutron scattering differs from, 
e.g., Ref. [15[ , where the pp-channel is considered the 
main contribution to the "41 meV resonance". Some fur- 
ther comparison of these viewpoints has been given in 
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FIG. 4. Wave-vector scan of x"(q> w ) (i n ar t>- units) 
around q = (71% 71") in the SC state at T — > . q x ,q y are in 
units of 2-7T =lr.l.u. The hole filling is x = 0.12 , energies are 
uj = 0.51J = ujq , 0.3J = 0.7ojo , 0.1J (from top to bottom). 
A quasi-particle damping F — 0.02J has been used. 

We now turn to the discussion of the magnetic re- 
sponse in wave-vector space. The resonance at 41 meV, 
as well as its relative in underdoped samples at energies 
luq < 41 meV , is characterized as a single (commensu- 
rate) peak at q = (7r,7r) = (1/2, l/2)r.l.u. ||-|f| . Re- 
cent INS experiments by Mook et.al. |L6| , performed at 
an energy loi = 24meV ~ 0.7loq on an YBa2Cu306.6 
sample with loq — 34meV gave evidence for an incom- 
mensurate response. It is dominated by four peaks, hor- 
izontally displaced from (ir, tt) . The maximum inten- 
sity is strongly reduced compared to uj = ujq . Fig. |] 
shows x"(l) w ) f° r x = 0-12 at T — * from the renormal- 
ized RPA Eq.([lj) with vertex corrections omitted. At 
the energy loq = 0.51J (top panel) the susceptibility 
is commensurate . When the energy is lowered to 
LJi — 0.7^0 = 0.3 J , the intensity drops dramatically, and 
four peaks appear at q = (tt ± 5, tt) and (tt, tt ± ft) , as is 
seen in the experiment [ fi6| . The amount 5 — 0.1 r.l.u. of 
the displacement fits well the observed value 0.105 r.l.u. 
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The incommensurate pattern may be characterized by 
the intensity ratio Ih /Id , which is the maximum inten- 
sity Ih found in a horizontal scan through (tt, 7t) , related 
to the maximum Id from a diagonal scan. The numeri- 
cal calculation yields Ih/ Id = 1-4; from the data given 
in Rcf. Jl6| we estimate an experimental value < 2.0 . 
When the energy is further reduced, Ih/Id changes con- 
tinuously to values < 1 , and for uj << ujq (bottom panel) 
four peaks appear at q = (tt ± 8', tt ± 8') [p6p|] . At these 
low energies no information on the q-space structure has 
yet been obtained from INS, due to the very dim inten- 
sity. 




FIG. 5. Left: contour plot of E(k) in the SC state for 
x — 0.12 . Shown is the upper right 1/4 Brillouin zone (BZ) 
with k x i y in units of 2-zr =1 r.l.u. The node E(k) = appears 
as a dot, encircled by contours E(k) = uj/2 for uj = 0.1J 
(most inner curve), 0.3J , 0.51J , 0.7 J . Right: contour for 
uj — 0.3 J and its by (tt, tt) shifted image in the full BZ . Arrows 
indicate the best nesting vector relative to (tt,tt) . 

An explanation for the incommensurate pattern at 
uj = uji can be found in the quasi-particle dispersion 
E(k) of the d-wave SC state. Consider again the ph- 
excitation energies f2(q, k) . At low energy uj — > only 
excitations between the nodes E(k) > , E(k + q) > 
contribute, with q = (tt ± 6', tt ± 5') near (n, tt) |gq] . For 
small uj = 0.1J this is still the dominant process, see 
Fig. H (left). The situation changes for uj = uji = 0.3 J. 
The contour E(k) = uJi/2 contains a flat piece, which 
allows for a nesting contribution from (nearly) degener- 
ate excitations E(k) sa E(k + q) = u>i/2. Following 
Ref. pjj the best nesting vector is a horizontal (or ver- 
tical) offset to (tt, tt), i.e., q = (tt ± 5, tt) and (tt, tt ± 8) . 
This is illustrated in Fig. |^ (right), the offset is given by 
5 = arcsin[(oj + fit)/ — t] . The above reasoning has been 
given without consideration of the renormalized RPA, 
Eq. ([!]) , since in fact x ~ X° f° r energies below ujq . On 
the other hand, this is not the case near the resonance 
energy ujq — 0.5 1J, where \ is determined by a pole in 
Eq. (^) , which produces the strong commensurate peak 
displayed in the top panel of Fig. || . 

In conclusion we find that the slave-boson mean- 
field theory, enhanced by the renormalized RPA, gives 
a reasonable explanation for the INS experiments on 
YBa2Cu306+j/ compounds. The energy and vanishing 
damping of the "41meV resonance" is associated with 



the threshold of ph-excitations at q = (tt, tt) . The evo- 
lution of the resonance with hole filling is accounted for, 
but the damping is quite sensitive to the bandstructure 
(i.e., t' 5^ 0). The same theory explains the incommen- 
surate structure at lower energies as a dynamic nesting 
effect, specific to the d-wave SC state. 
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